N PS  ARCHIVE 
1966 
FERRENTINO,  P. 


;;dffl 


BnnOOORi 
""MM 


WW, 


A  TWO-DIMENSIONAL  OMEGA  EQUATION  FOR 
THE  1000-700  MB  LAYER  WITH  DIABATIC  HEATING 


PETER  S,  FERRENTINO 


$3i 


ibSBffl 


ismmm 


i 


Bff» 


mi 


m 


1 


tfffifi 


matin 


m 


GEMS 


msM 


msam 


BuWMsHHbw 


■1 


i'v  ■"', ■  -  ''  I i 


m 


mw 


BOOK 

&$Be 

HI 

:  ''!.''                        ':  ■  -    '■  ■  '              '; 

WfflH 

Wlttl 

tnUifi 

|jj» 

■:i::;'.?;''v7i   ■■''''■'     ' 

rWP 

''''';■ 

'',"  ',* 

JK]B$||  .  \ 

B 

fi 

";'',■■'' 

V  'Nf'^U.:    ::         ..         .' 

mJtHV 

iJwfn 

H 

itnwi 
11 

V1' 

TJffiW 

■  'ii'i  :L     ''      !.■''■■■'■'...■ .    ' 

■■■ 


'■■'■''• '■;-::'  ;';',/ 


NHEK 


B5K 


ffissv 


■ 


SdCoct 


LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CALIF.  93940 


DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA   93943-5101 


This  document  has  been  approved  for  public 
release  and  rale;  its  distribution  is  unlimited. 


6 


&3 


NAVAL  POSTGRADUATE  SCHOOt 


A  TWO-DIMENSIONAL  OMEGA  EQUATION  FOR  THE 
1000-700  MB  LAYER  WITH  DIABATIC  HEATING 


by 

Peter  S.  Ferrentino 
Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment 
for  the  degree  of 

MASTER  OP  SCIENCE  IN  METEOROLOGY 

from  the 

UNITED  STATES  NAVAL  POSTGRADUATE  SCHOOL 

May  1966 


ABSTRACT 

A  two-dimensional  omega  equation  is  derived  by 
combination  of  the  vorticity  and  thermodynamic  equations. 
The  desired  omega  is  then  taken  to  be  the  logarithmic  average 
in  the  1000-700  nib  layer.   A  diabatic  term,  after  Laevastu, 
for  oceanic  areas  only  is  included  to  deduce  the  empirical 
temperature  and  vapor-pressure  changes  associated  with 
sensible  and  latent  heating  in  the  maritime  layers.   Over 
both  continental  and  oceanic  areas  a  frictional  vorticity 
sink  is  included  in  order  that  excessive  energy  cannot  be 
generated  over  the  ocean.   Among  other  novel  features  is  the 
use  of  the  Holl  static-stability  parameter  which  affords 
vertical  consistency  with   analyses  prepared  by  Fleet 
Numerical  Weather  Facility. 
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1.  Introduction 

Little  practical  us©  has  been  mad©  of  vertical  motions, 
CJ«  dp/</t,  in  short-range  forecasting  due  to  the  complicated 
and  lengthy  computations  required.   Vertical  motion  computa- 
tions  are  usually  the  by-product  of  multilevel  baroclinic 
models  and  single  values  must  be  extracted  only  after  going 
through  the  entire  three-dimensional  procedure.   In  this 
paper  a  two-dimensional  vertical  motion  equation  is  developed 
using  vertically  integrated  parameters  with  an  assumed  vertical 
motion  profile.   The  layer-mean  omega  is  especially  useful 
for  short-range  thickness  forecasts  and  the  1000-700  mb  layer 
has  been  selected  for  study.   Derivation  of  the CJ  -equation 
is  similar  to  that  by  Thompson  |_18J  except  a  diabatic  heating 
term  providing  a  mechanism  for  development  over  oceanic 
areas  has  been  included.   Due  to  the  empiricisms  employed 
the  non-elliptic  conditions  mentioned  by  Pedersen  [l6J  do 
not  arise. 

2.  The  Basic  Omega  Equation 

The  diagnostic  omega  equation  is  derived  by  standard 
means  from  a  vorticity  equation  and  a  thermodynamic  equation 
which  retains  the  diabatic  term.   Equation  (1)  is  the  time 
derivative  of  the  first  law  of  thermodynamics  in  (x,  y,  p,  t) 
coordinates. 


ii 


Substituting,  T  =-|  $?     ,    which  comes  from  the  hydrostatic 
assumption  and  the  ideal  gas  law;  then  dividing  by,  -g/R, 
yields  equation  (2). 


•pfitlno-  on  enuation  ( ?)    with  the  Laolacian  gives 


(2) 


Operating  on  equation  (2)  with  the  Laplacian  gives  the  final 
form  of  the  thermodynamic  equation  as  shown  by  equation  (3). 

Here,  O^s— Rp3*^,  is  the  Holl  stability  parameter  to  be 
further  discussed  in  section  (i|). 

The  vorticity  equation  in  pressure  coordinates  is  shown 
by  equation  (1|) . 

As  presented  by  Thompson  llSjand  numerous  other  writers 

the  last  term  of  the  left  side  of  (I4)  is  approximately  equal 

to  that  on  the  right  side,  and  the  two  terms  are  henceforth 

deleted.   Then,  a  useful  form  of  the  vorticity  equation  is 

arrived  at  (equation  5)  by  making  the  geostrophic  assumption 

for  vorticity  and  for  velocity,  C  -  JLV*  >    W,*$.WZf   and  by 

0  -   f  *  f 

taking  the  logarithmic  pressure  derivative. 
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Here,  7*5^1  >  is  the  absolute  vorticity.   Next  subtract 
(5)  from  (3),  and  the  result  is  the  omega  equation: 


-T8^'^)) 


(6) 


3.   Vertical  Distribution  of  Pressure  and  Omega 

Instead  of  the  distribution  of  pressure  indexes   normally 
considered,  that  employed  here  is  given  by, 


where  V)    is  an  arbitrary  pressure  level.   The  1000-500  mb 
layer  is  divided  into  sub-layers  as  shown  by  Figure  (1). 
The  resultant  pressures,  while  essentially  logarithmic,  bear 
values  nearly  identical  to  those  of  the  mandatory  levels. 
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Figure  1.   Pressure  Distribution  in  the  Vertical 
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The  values  of  OO-  dp/dt,    the  "vertical  velocity  in  pres- 
sure coordinates  ,   La  assumed  to  vary  parabolically  in  the 
vertical  with  a  profile  defined  by  equation  (7)  which 
extends  nearly  to  the  5>00  mb  level. 


u>=A(j<ALt+bfA2.\ 


(7) 


For  boundary  conditions  it  is  assumed  that^OOat  p=p0 
(the  kinematic  boundary  condition  at  the  assumed  surface  of 
the  earth,  p0  ).  Also  for  convenience  it  is  assumed  that  at  p 
C/C^a  0  ,    i.e.  that  the  3>00-mb  divergence  is  zero.  Then, 
upon  differentiation  of  equation  (7)  with  the  500-mb  diver- 
gence taken  to  be  zero,  equation  (8)  results.   It  should  be 
noted  that  if  the  value  ( ^%yp)«0  is  assumed  zero  at  p=p3 
the  value  of  B  is  altered  by  only  5$. 


or. 


Therefore,  B=1.3863A,  and  the  omega  profile  in  terms  of  A  is 
given  by  equation  (9). 


cb-mA 


LA  +  i.m3^ 


(9) 


Since,  CO 0  =o,  the  ratio  between  60,,  and  ^)^,  defines  the 
omega  profile  (equation  10)  for  the  0-2  layer  (see  Figure  2) 


Ik 


CO, 


^AAz(- £/**+/.  3863) 


(10) 
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Figure  2.   The  Omega  Profile 


For  the  pressure-range  considered  (see  Figure  2),    the  layer- 
logarithmic-mean,  UJ  =0.950dO|,so  that  U>, ,  will  be  used  f  or  U)  . 
Note  finally,  the  profile  assumed  here  applies  only  to  the 
large-scale  adiabatic,  frictionlcss  component  of  the  vertical 
velocity  (that  is  terrain  irregularities  are  considered 
absent ) . 

l\.      Vertical  Integration  of  the  Omega  Equation 

The  Holl  stability  parameter  Lllj  is  used  since  it  is  one 
of  the  vertical— consistency  requirements  used  by  FNWF  (Fleet 
Numerical  Weather  Facility)  and  this  study  uses  FNWF  processed 

data.   The  Holl  parameter  is  a  modification  of  the  standard 
stability  parameter  used  here,  as  shown  by  equations  (11),  (12), 
and  (13). 
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&dx>         r\Qp?      dp  J 


c?f      ap 


(id 


Letting,  RT=  -  <3j  ^|L    gives, 

aS£Tp 


(12) 


Equation  (12)  may  now  be  finite  differenced  over  the  0-2  layer. 


<r« 


§-(**-*•) -(To-Tl) 


(13) 


Equation  (13)  represents  the  Holl  stability  parameter  in 

finite-difference  form.   In  effect, Oh  ,  by  (13)  gives  the 

i, 
logarithmic-pressure  \ 'average'    over  the  0-2  layer.   This 

- 

use  of  O^  is  compatible  with  the  assumed  logarithmic  pressure 
distribution  and  will  be  used  throughout  the  analysis  which 
follows. 

The  omega  equation  (6)  will  now  be  integrated  term  by 
term  over  the  isobaric  layer  (Po,  ?z) •  Prom  (6),  one  may 
rewrite  the  QJ -equation  as: 


,/(D)   /     v\        (E)  X' 

whose  parts  (A),...,,  (E)  will  now  be  discussed  individually. 
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(A)  The  first  term  is   approximated   by  substituting 
layer  mean  values:   ^^60,)  ,  60,  =  U) 

(B)  With  the  mathematical  relation,  p*  £> <~Q  -  &  Q        h<*>  r 
the  second  term  becomes,     .  J    \  X*  ■ 

whose  center-finite-differenced  form  over  the  0-2  layer 


becomes: 


(ii+) 


2 

Then,  using  the  specified  vertical  motion  profile  the  rela- 
tion^2.=  1.711+3^\  is  substituted.  However,  it  is  desired 
to  introduce  lower  boundary  effects  due  to  friction  and 
terrain  by  letting  &Ja»£0|.o  ,  the  U)uo  term  then  being  pro- 
determined  and  placed  on  the  forcing  function  side  of  the  CO- 
equation.  Expression  (15)  below  is  the  final  form  of  (11+) 
including  effects  of  the  lower  boundary. 

Term  (B)  i  +%  (     ^M2JSi     |  A?/V3*>,\ 

j.fW/_ai ^o  ^    (15) 

The  lower  boundary,  LOUO  ,  will  be  discussed  in  section  (5). 

(C)   Term  (C)  in  the  geostrophic-diagnostic  model  repre- 
sents geostrophic  advection  of  absolute  vorticity.   Averaging 
in  the  vertical  with  respect  to  the  logarithm  of  pressure  leads 
to  equation  ( 16) : 
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j-  MjM)i^pi  m  (ajU^Iky       (i6) 

The  derivation  of  equation  (16)  makes  use  of  the  mean  value 
theorem,  so  that  p  =  p-  is  represented  by  p  and  may  be  re- 
moved from  within  the  integral.  Equation  (16)  then  becomes, 

M.l(*,ty  .iil^^-Jil^l)  d7) 


Lty. 


(D)   Term  (D)  involves  the  geostrophic  advection  of 
thickness  within  the  logarthmic  pressure  integral.   Its 


value  is  well  approximated  by  equation  (18): 


ftv##'ih»' 


the   right  side  of  which  gives  equation  (19). 

h£vzJ(z>  ,a*)     « Ii£  ^J(e<,  zz-so)  d9) 

Equation  (19)  shows  that  the  geostrophic  advecting  wind  may  bo 
taken  arbitrarily  as  the  level  (1)  wind. 

(E)   The  diabatic  heating  term  will  be  vertically 
integrated  in  section  (6). 

The  integrated  form  of  the  resulting  omega  equation,  subject 
to  the  lower  and  upper  boundary  conditions  (atpfcandp^_) 
becomes, 
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+t  = 


Equation  (20)  is  a  Helmholtz-type  equation  in  the  variable, 
5L60,  .   This  grouping,  0*^60,  ,  precludes  making  the  usual 
simplifying  approximation  of  most  three-dimensional  models. 
(see  for  example,  Haltiner  et  al  9  [I. 

F2^  to)  -  <rv*a)  -i-  6o\v V  +  z«/(r-  W6o  &  cnvio 

The  final  solution  of  equation  (20)  is  to  be  divided  by  O^ 
leaving  U)t=  ^O  ,  where ^  is  the  mean  vertical  motion  which 
approximates  the  8£0  mb  vertical  motion. 

5>.   The  Lower  Boundary  Condition 

Vertical  motion  at  the  lower  boundary  (Muo,    is  the  sum  of 
a  terrain  and  a  friction  term,  Lduv  a  ^T-\-  ^F  • 
Here,  CO-f   ,  the  terrain-effected  vertical  motion  may  be  de- 
scribed by  equation  (21)  from  Berkofski  and  Bertoni  [  2  [  . 

Cor  «  -ftft/r  •  IP^t  (21) 

Terrain  height,^,  is  a  smoothed  field  used  operationally 
by  FNWF.  The  wind  velocity  at  terrain  level,  IV-p  ,  will  be 
arrived  at  by  an  objective  procedure  to  be  described  below. 
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The  friction  term  L0?  ,  as  shown  by  equation  (22)  is  a 
simplified  version  of  Gressman's  formula  used  by  Haltiner 
et  al|"9]  . 


CO?  =-f>T3L  Ct>Vt$t 


(22) 


Q>  is  the  geostrophic  drag  coefficient  derived  by  Cressman  j  7  [> 
the  field  of  which  is  in  regular  use  at  PNWF,   The  field  of 
C~£  is  a  set  of  constants,  one  for  eachAj]  grid-point  in  the 
Northern  Hemisphere. 

Since  the  J"   subscripts  infer  application  at  terrain  height, 
values  are  used  that  approximate  the  particular  terrain  heights 
involved.  This  is  done  by  dividing  the  0-2  layer  into  three 
terrain-height  dependent  cases. 


_  A  _700_  m 


1  J00_m_ 


?z  Case   3:    2200   m  £  2-r 

-2200  m  use   Pz  values 

Case   2:      900  m^Z/2200  m 

900  m  use   Pi    values 

Case  1:      sfc^&r^OO  m 
use   P0  values 


p 


Figure  3.   Three  Cases  of  the  Lower  Boundary 


Division  of  the  0-2  layer  as  shown  by  figure  (3)  is 
dependent  upon  the  gradient  level  above  terrain,  since  \\/T    is 
approximated  by  some  geostrophic  wind  throughout.   As  a  basis 
for  the  selection  of  the  particular  geostrophic  wind, \\Aj  ,  note 
that  the  gradient  level  in  a  neutral  atmosphere  occurs  at 

2| 


approximately  700  m  above  terrain  level  for  a  wide  range  of 
surface  roughness,  as  shown  by  Blackadar J  3   .   Therefore, 
vertical  velocity  at  the  lower  boundary  is  case  (1)  if  the 
terrain  height  is  900 m or  lower.   This  values  corresponds  to 
a  "standard  atmosphere"  height  which  is  700m  or  more  below 
the  (1)  level,  and  (0)  level  parameters  are  then  used  to 
compute  the  lower  boundary.   Case  (2)  is  to  be  used  if  the 
terrain  height  j?T  lies  within  a  gradient-level  range  of  pf 
but  not  P^  .   Case  (3)  occurs  for  all  terrain  heights  within 
700  m  of  pt  , 

Equation  (23)  is  the  terrain-induced  vertical  velocity  at 
the  lower  boundary  after  application  of  the  geostrophic 
assumption. 


> 


(coj)^  =  _  a,  £.  T(b* ,  2T)  -^  £  4 1  'i<*  t>»  v\  <  23 

Here,  the  V\      subscript  refers  to  the  levels  0,  1,  and  2. 
Densities  are  standard  atmosphere  values  for  the  respective 
levels, 

6.   The  Diabatic  Heating  Term 

Diabatic  heating  in  this  model  results  from  an  exchange 
of  heat  between  the  sea  surface  and  the  atmosphere.   The  ex- 
change is  divided  into  three  cases;   neutral,  convective,  and 
inversion,  the  division  being  dependent  upon  the  sea-air 

temperature  difference,  T  ~  T  . 

*     w    a 

The  neutral  case  (0^  C  £  T.  -TA<3°C)  results  in  no  net 
transport  of  heat  across  the  air- ocean  interface  and  the 
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diabatic  heating  term  is  arbitrarily  taken  as  zero  when  this 
condition  prevails. 

^fe_^g_gYe-C-^Ye_Jlig£  (  i^v"ta-  3°C)  results  in  a  net  heating 
of  the  atmosphere.   The  oonvective  atmosphere  is  assumed  to 
consist  of  a  dry  adiabatic  lapse  rate  to  the  condensation 
level  and  a  moist  adiabatic  lapse  rate  above.  This  particular 
model  of  the  convective  case  is  discussed  by  Burke  [5 J  in  his 
paper  on  the  transformation  of  cP  to  mP  air.  Heat  gain  by  the 

atmosphere  is  manifested  by  a  layer-mean  temperature  gain 
induced  by  a  one-hour  trajectory  of  surface  air  over  a  warmer 
sea-surface,  and  the  layer  temperature  change  is  shown  by 
figure  (1|). 
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w                                                        "■^^"«                                       \'  ■ 
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\&fo\ 

A U 
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Figure  k»      Model  of  the  Convective-Atmosphere  Modification 


Heating  is  composed  of:   (1)  sensible  heat  from  the  sur- 
face to  the  condensation  level  and  (2)  latent  heating  above 
this  level.  This  result  follows  by  an  analysis  of  the  Burke 
model  and  its  implications  reguarding  heat  transport.   Martin 
I  lljj  has  shown  that  the  large-scale  transports  of  heat  have 
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already  been  included  In  those  terras  of  equation  (6)  in 
which  Q  does  not  specifically  appear.   Thus  Qi/Cp  ,  for  use 
in  equation  (6)  is  given  by  the  local  temperature  change. 

Equation  {2$)    is  the  diabatic  heating  ten  (■£)  referred 


to    '  hat  ion    (6).. 


d?w^A  --Rpv* 


ATa\  _  /ATwa\ 


At  / 


(25) 


Integrating  logarithmically  in  the  vertical  yields  the 
final  form  of  the  diabatic  term  as  shown  by  equation  (26). 


(26) 


In  each  of  equations  2fy,  25,  and  26,  the  operator-^  represents 
the  finite-difference  version  of  the  local  time  derivitive, 
while,  P^  ,  represents  a  convective-condensation  level 
pressure. 

According  to  the  model  depicted  in  figure  (1|),  the  modi- 
fication is  largely  determined  by  the  static  stability  which 
in  turn  largely  depends  on  OJ  .  As  an  assumed  infinite  source 
of  heat,  the  ocean  will  first  modify  a  thin  surface  layer  of 

air  and  then,  due  to  convective  activity,  will  distribute  the 
modification  (temperature  change)  throughout  the  0-2  layer. 
An  overestimation  of  heat  exchange  should  be  expected  with 
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this  particular  process  since  convection  is  not  instantaneous 
This  overestimation  may  be  partially  balanced  by  heat  lost 
from  the  top  of  the  0-2  layer.  However,  the  non-inclusion  of 
a  moisture-continuity  equation  may  be  more  serious  in  regions 
of  upper  ridges. 

The  inversion  case  (T^  -  TA<  0)  or  c_cjolin£  o^sej.  where 
heat  is  lost  from  the  atmosphere  to  the  ocean  is  defined 
by  an  empiricism  from  Laevastu  [ 12J. 

Qc  =  3,0  Vic   (Tk;  -Ta)  gm-cal/cmx-2Lthrs         (27) 

Equation  (27)  represents  the  heat  flux  across  the  10-meter 
level.   This  is  an  empiricism  requiring  the  10-meter  wind, 
V|o,  which  will  be  approximated  by  making  a  frictional 
correction  to  the  1000  mb  geostrophic  wind  as  shown  by 
equation  (28) : 

V)a«  a.sJL  IkxvDo  (28) 

f 

The  factor  0.5  is  due  to  surface  frictional  effects  and  will 
be  discussed  in  section  (7). 

Before  integrating  equation  (27)  over  the  0-2  layer 
the  initial  distribution  within  the  layer  must  be  determined. 
Normally  cooling  effects  due  to  inversion  conditions  are  con- 
fined to  a  surface  layer  only  several  hundred  meters  thick, 
however,  since  the  layer  mean  cooling  effect  is  desired  over 
the  entire  0-2  layer  it  will  be  assumed  that  cooling  is  dis- 
tributed throughout  the  0-2  layer.   Thus,  dividing  equation 

2h 


(27)  by  the  mass  of  a  column  of  air  extending  from  the  surface 
pressure  PA  tc  Pa,  with  total  mass,M  =  (PA  -  PjJ/g,  gives  the 
mean  temperature  change  in  a  vertical  column  extending  from 
PA  to  Pt  as, 

A.  -.      ®*3L  (29) 

Finally  we  have,  ?a, 

_  RaV*£  =  -  El^lF^    ^,  -  £*-  V"q       (30) 

which  is  the  final  result  for  the  inversion  case  assuming 
neglible  heat  release  from  fog-droplet  condensation, 

7.   State-Change  Parameters  for  the  Convective  Case 

Air-temperature  change  and  wet-bulb  temperature  change 
equations  must  be  derived  for  use  in  equation  (26).   The 
state-change  equations  (31)  and  (32)  were  initially  developed 
by  Mosby  [l5L  Amot  IJ,  and  more  recently  by  Boyumjij  J. 
However,  in  a  FNWF  study  by  Carstensen  and  Laevastu[6 J  the 
following  best-fit  equations  were  found  to  give  hourly  changes 
with  a  high  degree  of  skill: 

£Ja  =.aiz.(%j-T^-o.io  1'vio'B^  °cAr       (31) 
At 

£*  *  o./5{eui-eA)-  0J6  -  W/VFe*    mt/hr        (32) 
it 

Equation  (3D  is  substituted  directly  into  the  sensible- 
heat  term  of  equation  (26)  while  the  wet-bulb  temperature 
change  contributes  directly  to  the  latent-heat  term.   In 


25 


computing  AJ^l  it  is  necessary  to  Include  the  small  Increment 

At 
of  temperature  found  by  proceeding  along  a  mixing  ratio  isopleth 

from  P^,   to  ?c^   then  down  a  moist  adiabat  to  Pc,  (see  figure 

5).   Equation  (33)  shows  this  relationship. 

ATwa  =  &Ti>-t(f*-Xn)AZ}*-  ,  &&*  J.406-2L   (33) 

Here  A^pc  is  the  one-hour  height  change  of  the  condensation 
heights  determined  by  the  height  of  intersection  of  the  dry 
adiabat  and  the  mixing  ratio  isopleth  based  upon  the  tempera- 
ture -  dew-point  spread  at  the  surface. 


A2pc  =  2>ci  -2 


Figure  5.  ^TL,o  as   a  Function  of  /IVTk 


(31+) 


Substituting  equation  (3k)    into  (33)  gives  the  local  rate  of 
change  of  the  wet-bulb  temperature  as  a  function  of  dew-point 

and  a i r - 1 o nrnorn tnro  ohnnro, 


A"U  =  AT*  ,  //d-W)/aTa  _  A/b 
At       At       ^t>-tfd  /I  At       At 


(35) 
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The  dew-point  temperature  change  /$i  is  found  by  taking 
the  time  differential  of  equation  (36)  below  10 1, 

TA-Th  =  a^L(es-ga)  06) 

Here,  es    is  the  saturated  vapor  pressure  and  e^  the  observed 
value,  both  at  10-meters.   The  finite-difference  form  of  the 
time  differential  of  dew-point  resulting  from  equation  (36) 
then  follows: 

^j^^Alji  .  £x^zLAe*      a^4^^LAeg         (37) 
At.      At-       4>$    At  <t?>    &< 

The  vapor  pressure  change,  &&A     t    is  given  by  equation  (32). 
The  time  rate  of  change  of  saturated  vaoor  pressure,  ^5  , 
is  calculated  using  the  Cla&sius-  Ol'afjfyron  equation: 

^^  KTa^    At  (38) 


fir    6JO?<fXp 


$<//*,£    SV/?,o 


(39) 


293.  \(*  7> 

Equation  (39)  is  an  integrated  form  of  the  Clausius  -  Clapyron 
equation. 

There  are  still  some  unsolved  parameters  in  the  system. 
The  moist  adiabatic  lapse  rate,  ^ ,  cannot  be  considered  a 
constant  with  respect  to  pressure  even  in  the  limited  opera- 
tional range  of  this  study.   At  the  condensation  level,  Pc  , 
the  moist  adiabat  is  calculated  from  the  following  |"lo"j  : 
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£>,4>ZZ  L<2 


(1(0) 


All  subscripts  G  refer  to  the  condensation  level. 

Condensation  parameters  of  temperature,  Tc,  and  pressure, 
Pc  ,  are  readily  calculable  from  equation  (l[l)    and  (lj.2)  after 
Edson  j^8]: 


to- Si 


Pc 


=K*)% 


(III) 


(1*2) 


For  equation  (1|0),  the  vapor  pressure  at  the  condensa- 
tion level,  ec,  must  be  calculated.   Equation  (k3) »    the  inte- 
grated form  of  the  Clausius  -Clapeyron  equation  may  be  used 
since  air  at  the  condensation  level  is  saturated,  and  Pc  and 
Tt  are  known  from  (1|1)  and  (1|2), 


£<.  J&.lofetp 


272.1(0 


T, 


(1*3) 


In  like  manner  the  vapor  pressure  of  the  sea  surface,  e  , 
(from  equation  32)  is  calculated  since  the  sea  surface  is 
saturated  (equation  l\)\) ,      Also,  a  correction  factor  of 
0.98  is  required  for  salinity  effects. 


UJ 


=  (b.«?  )(4/o?)s><f> 


Z93./4 
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Finally,  the  advecting  wind  at  lO-meters,^^,  must  be 
derived  for  use  in  equations (31)  and  (32).   The  geostrophio 
wind  at  the  surface  is  taken  as  equal  to  the  1000  mb  geo- 
strophio wind.  However,  due  to  frictional  interaction  with 
the  ocean  surface,  the  geostrophic  wind  must  be  reduced  in 
magnitude  and  rotated  to  the  left,  to  yield  the  advecting 
frictional  wind,!J/,0.  For  the  convective  case  the  observed 
cross-isobar  angle  at  mid-latitudes  is  approximately  lf>  L10!* 

Reduction  of  the  10-raeter  wind  speed  as  expressed  by 
the  ratio,  V/V*0  ,  where  V^  is  the  10-meter  geostrophic  wind 
speed,  may  be  obtained  from  a  study  by  Lettau  |  13 "I.   By 
assuming  a  roughness  length  of  £0  =  0.1  cm  for  oceanic  areas, 
the  10-meter  geostrophic  wind  ratio  for  neutral  conditions 
is,  V/V*0  =  0,6(i   Lottau  then  allows  for  stability  criteria, 
and  for  the  convective  case  a  value  of  V/V«0  =  0.7^' is  arrived 
at.   Similarly,  V/V*g  =  0.J>  for  the  inversion  case. 

Advection  by  the  frictional  wind  is  derived  by  first 
decomposing  the  vector  advection  as  shown  by  equation  (1|5>): 

Then  multiplying  by  O.1^  and  rotating  the  field  through  the 
angle  K  =  l£  by  a  rotational  change  of  coordinates  leaves 
equation  ( Ij.6 )  * 
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Assuming  the  geostrophic  conditions,  Ufo*  —£  S^>  ^--zr^ 
equation  (i|6)  may  be  written  in  ita  final  form  aa  given  by 
equation  (i|7). 


(1*7) 


Similarly,  equation  (!j8)  is  derived  for  the  adveotion  of 
vapor  pressure  by  tie  10-meter  wind. 

V>*'We*  =  0.1  l_<u*ttaT  (*■<>£  a) 

f 

8.  Numerical  Procedurea 

All  finite-difference  operators  used  are  standard  cent- 
ered-differencea  of  the  type, 

V^A  =  2aIVM  J(A,b)  =  JkL.  J(atb) 

with  five-point  grida  having  a  mesh  distance,  d  =  3^1  km 
at  60   latitude.   !Iere,M  is  the  map  factor  for  polar  stereo- 
graph projections  and  d  is  the  grid  spacing. 

The  variables  needed  for  input  into  the  diagnostic  omega 
equation  x^ere  scaled  as  follows: 
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A  A         l! 

m      =>Vt  .2                                  p  =p.  .2     rab 

)9       *y)   .2"  sec"1                     €  =C  .2**  mb 

^O      =60.2  mb/sec                 CJJ  =0J.2     crnVsec* 


-A        ?        /  —         4        *° 


cm 


^Jto  =  ^.2      mb/sec  ^T  =  :2t*2 

}>       =b  .2/?  cm  ^    =<$    .2"    gm-cm7sec* 

T      =|     .2     C  -T    =     l.lj  £6142x10"    sin^sec-' 

The  following  physical  constants  were  also  needed  in  the 
computations  of  Oj  .   All  values  of  those  constants  have 
been  expressed  in  cgs  units. 


d 

= 

38,1000,000   cm 

* 

— 

980.0   cm/sec1 

R 

s 

.0287  X  107    cmVsec* 

-   °K 

K 

as 

0.2+61  X   107    cmx/secl 

-    *K 

c* 

= 

1.003  X   107    crnVsec1- 

-  *K 

L 

= 

2500  X   107      crnVsec* 

oT/- 

-       ft. 

*d 

= 

0.9771  X   10_t/  °C/crn 

rfc 

= 

0.1600  X   10        C/cm 

/* 

= 

1.213  X   10"3    gm/cm3 

c> 

ss 

1.0££  X   10" 3    gm/cm* 

r* 

= 

0.919  X   10"*    gm/cm3 

The  scaled  Helmholtz  equation,  which  follows  from 
equation  (20),  expressed  in  units  of,  mb  -  cm2/sec3  ,  is1 
shown  on  next  pace. 
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vt^A)-  H.^osilhi^cb)'^  \ 


<rH» 


x&^-kX,/^ 


.1-/ 


A 
"2  ♦ 


'2  -p^^.Z 


-42 


049) 


All  height  values  have  been  replaced  by  D-values,  the 
deviation  of  the  height  field  from  a  standard  atmosphere 
height. 

The  scaled  stability  parameter  is  given  by  equation 
(5>0)  and  is  in  units  of  cmI/sec^ 


■  Z'taooo 
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^Pt^xr  (5o 


For  computation  of  the  stability  parameter  D  is  re- 
scaled  to  D  =  D   .2   to  prevent  overflow  when  adding  the 
standard  height, 2>6~^/000  ■  290,000  cm.   Stability  calcu- 
lations are  based  on  1000-700  mb  differences  and  throughout 
the  derivation  all  gradients  and  differences  at  the  P,  and 
Pz  levels  are  assumed  equal  to  those  at  8£0  mb  and  700  mb 
respectively. 

To  prevent  computational  instability  during  the  numerical 
Helmholtz  solution  a  minimum  value  of  stability  is  required. 
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This  value  corresponds  to  a  maximum  lapse  of  ~%X(\ : 

(fo-7.)  *"$*£  [(M'}?"  +  21°°°  c  " 

Equation  (51)  is  the  scaled  lower  boundary  condition 
in  units  of  mb/sec. 


(5D 


Here,  the  subscript  Y\  refers  to  the  terrain-height  dependent 
cases,  1,  2,  and  3  explained  in  section  (5).  The  factor  10"" 
converts  dynes/cm2-  to  millibars. 

Equations  for  the  diabatic  term  are  computed  in  terms 
of  heating  rates,  Q/Cp  ,  which  are  assigned  to  grid  points 
according  to  the  existing  grid  condition:   convective, 
inversion,  or  neutral. 

Equation  (52)  shows  the  scaled  heating  rate  for  the  convective 


i* 


(52) 


case. 

A 

A  - 

The  factor,  3600,  converts  hours  to  second  and  heating  rate 
is  now  in  units  of  mb  -  aC/sec. 

Equation  (53)  shows  the  scaled  heating  rate  for  the 
inversion  case  (a  cooling  process)  also  in  mb  -  °C/sec. 
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J  A)  =  4A  fa^^^^?^^-^  ,53 


Cp 


f2fl( 


ICC 


f 


The  factor  (i|.l81|  X  10)  converts  gm-cal/sec  to  dyne-cm/sec; 
K  =  2\\   X  3600  which  converts  2l|-hours  to  seconds;  0.03  re- 
places 3.0,  and  permits  velocity  computation  to  be  made  in 
cm/sec. 

Land  area  grid  points  are  masked  and  receive  heating 
rate  values  of  zero  as  do  oceanic  grid  points  which  satisfy 
the  neutral  conditions.   The  remaining  oceanic  grid  points 
receive  either  heating  or  cooling  values  depending  on  the 
prevailing  condition  at  the  point.   The  diabatic  term  is  then 
computed  and  the  scaled  form  in  units  of  mb-cm  /sec3  is  shown 
by  equation  (ft) . 


(ft) 


Here,  P0  approximates  P  for  ease  of  computation. 

Other  equations  involved  in  the  computation  of  the  dia- 
batic term  are  scaled  as  follows: 


(1)   The  air  temperature  change  is  units  of  °C/hr. 

At  '  \   / 


-q 


-0.7M1 


^d1 


CX>StS^(4aJ^^^\Sm(^J^  +A^tf*\ 


(55) 


»* 
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(2)  The  dew-point  temperature  change  in  units  of  °C/hr. 


—  0.7A2H- 


d(&i%"T(ia ,  Ga)  -  siw«*(4$»A$+ 


+  ^l 


A^)] 


l^ 


Co.&Z2-L?'£.s)  AJfi 


Ai 


.X 


(56) 


(3)      The  wet-bulb   temperature   change    in  units   of  *C/hr. 


where, 


(58) 


^-0.497*  ^±!±^i!±^W 

(l±)     The  condensation  level  temperature  in  units  of  C. 

(5)  The  dew-point  temperature  in  units  of  C. 

(6)  The  condensation  level  pressure  in  units  of 


(59) 


millibars,     Pc  =& 


"& 


and 


(60) 


Hero  Tai    is   computed   using  values   of   Tft   and  T*    in  enuation  58. 


R*  =  Pa 
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Here  T~  is  calculated  using  values  of  TA|  +AT*/At  and 
ir+A^/At  in  equation?*  ( £6 ) ,  (59.),  (55),  and  (56). 

9.  The  Computer  Program 

The  Control  Data  Corporation  l60l|  digital  computer 
was  used  in  this  study.   It  has  a  core  capacity  of  32,768 
words  of  I48  bits  each.   An  operational  field  of  1,977 
grid  points  forming  a  5l  X  \\7   octagon  inscribed  within  the 
9°N  latitude  circle  was  employed.   The  grid-mesh  is  38I  km 
true  at  60°N  latitude. 

Boundary  conditions  around  the  octagonal  grid  were 
determined  by  the  standard  subroutines  used.   Laplacian, 
Jacobian,  and  all  other  five-point  center-difference 
operations  set  the  grid  boundary  to  zero.   The  Helmholtz 
relaxation  operation  set  the  edge  and  the  next  interior 
border  point  to  zero. 

Equation  {\\%   was  solved  by  a  two-dimensional  Liebmann 
relaxation  technique  wherein  the^l  +l)-iterate  for  any  point 


A 


is  givenby. 

A  *+l   /  " 


V^,;-(^b).i-C,-,- 


2-lB;,-  +  / 

Here,  n    is  the  over-relaxation  coefficient  and  the  residual 
at  any  step,  Rn  ,    is, 


4 


V^.i-(AB)ij-£'i 


-1  n 


-i 


t>r,  + 1 
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The  over-relaxation  coefficient  used  was,  }\  =  1.1*H[, 
which  allowed  convergence  in  approximately  30  scans  over 
the  1977  point  grid  using  an  initial  guess  field  of  zero.  The 

Convergence  criterion  is  that  the  iteration  ceases  when  € 

(n)  (v\±t\  (*)  z..   a 

defined  by  %    sA  -A  falls  below  3000  mb- cry  sec  . 

This  value  corresponds  to  a  vertical  velocity  of  0.3  X  10 

mb/sec  for  a  stability  of  83  X  10  cm  /sec   (the  standard 

atmosphere  stability). 

Total  computation  time  was  approximately  1  minute  and 
30  seconds  with  23  seconds  required  for  the  heating  term 
and  26  seconds  for  the  boundary  condition. 

Due  to  the  abrupt  cut-off  criteria  used  for  delineating 
the  three  cases  of  diabatic  heating  (inversion,  convective, 
neutral)  the  final  Q/e0    field  was  smoothed  with  a  five-point 
smoother  of  the  form, 

A  =  A  +  KSPA 

The  smoothing  coefficient  k  is  a  constant  value  of  (1/8) 
and  the  field  was  smoothed  twice.   The  smoothing  operation 
removed  small  scale  irregularities  caused  by  the  cut-off 
criteria  of  the  heating  term  but  also  reduced  peak  values. 

A  standard  FITWF  filtering  process  followed  all  operations 
involving  the  Laplacian.  Tills  process  removes  all  small 
scale  features  with  wave  numbers  greater  than  15  at  latitude 
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10.  Results  and  Conclusions 

A  series  of  four  successive  12-hourly  data-sets  be- 
ginning with  the  00Z  maps  of  28  April  1966  and  ending  with 
12Z,  29  April  1966  were  used.   Since  the  only  diabatic 

heating  mechanism  introduced  into  this  study  was  that  aris- 
ing by  "conduction"  from  the  underlying  oceanic  surface, 
the  results  have  been  depicted  only  over  the  North  Atlantic 
Ocean.   This  region  has  a  large  density  of  reporting  ships 
and  thus  the  reason  for  its  selection. 

The  computational  omegas  and  their  associated  parameters 
for  00Z  28  April  1966  are  contained  in  figures  (6)  through 
(11)  which  are  appended.   Of  these,  not  only  the  diabatic 
omegas  are  shown,  but  also,  the  additive  effect  of  the  dia- 
batic influence  as  contrasted  with  adiabatic  computations. 
In  figures  (12)  through  (1?)  only  the  final-product  omegas 
of  this  study  are  shown  together  with  the  corresponding 
FNWF  analyzed  8^0  mb  D-fields  to  servo  as  i  dent  if  dor's",: 
that  is,  to  indicate  qualitative  coherency  between  the  verti- 
cal motions  and  the  associated  motion  systems. 

In  the  sequence  of  figures  (6)  through  (11),  it  is  of 
interest  to  note  the  pattern  similarity  between  comparative 
adiabatic  gSJgJiAfi&l  motion  computations  rthose.  produced  here, 
and,  for  the  same  timrs,  those j  produced  tyy  F1F./F,  which  are  based 

upon  the  procedure  described  by  Haltiner  et  alf"9l.   The 
magnitude  and  position  of  the  updraf t-centers  show  a  strong 
similarity  with  those  of  FITV/F,  however,  the  model  described 
by  this  paper  shows  larger  downdraft  areas  extending  east- 

38 


ward  from  Newfoundland.   This  apparent  discrepancy  may  be 
attributed  to  the  two  differing  treatments  in  ^LO(figures 
18  and  19).  The^c^at  the  continent-edge  essentially  becomes 
the  boundary  condition  for  the  oceanic  computations. 

On  the  other  hand,  however,  one  of  the  major  objectives 
of  this  study  was  to  test  a  simplified  ^>l»  which  uses  only 
parameters  pertaining  to  the  standard  levels  of  analysis 
(equation  23).   In  this  connection,  the  greatest  time- 
consuming  aspect  of  the  FNWF  <-o  -computation  is  that  of 
deriving  U)uo,    and  associated  parameters  at  terrain  height  re- 
quiring pressure  extrapolation.   The  general  similarity  in 
the  fields  of  ^(.©by  the  two  operational  procedures  which 
have  been  discussed  is  evident  by  an  inspection  of  figures 
(18)  and  (19),  and  this  can  perhaps  be  a  justification  for 
continuing  the  simpler  ^©computations. 

The  diabatic  effects,  which  were  obtained  for  the  00Z 
April  28  synoptic  time,  depict  values  of  Q/up,  so  that  the 
large  positive  center  southeast  of  Newfoundland  is  realistic 
with  regard  to  the  northerly  flow  passing  over  the  Gulf 
Stream.  The  maximum  heating  value  for  the  four  map  periods 
over  the  Gulf  Strewn  Is  h/^v   =  O.O67  mb  -  °C/sec,  which 
corresponds  to  a  layer  mean  heating  rate  of  0.3  °C/hr  for  the 
1000-700  mb  layer.  The  equivalent  thickness  chaiige  for  this 
layer  is  approximately  3  meters/hr. 

For  comparison,  Petterssen»s  values  \yf\   for  the 
1000-500  mb  thickness  change  during  a  similar  synoptic 
situation  in  late  March  show  a  maximum  value  of  6  to  8 
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meters/hr  in  the  vicinity  of  the  Gulf  Stream.  These  values 
for  sensible  and  latent  heating  were  computed  using  flux 
calculations  at  the  surface  based  on  empiricisms  similar 
to  those  of  Laevastu.   It  should  be  noted  that  heating  rates 
for  this  time  period  in  late  March  are  indeed  greater  than 
late  April  values. 

Finally,  with  regard  to  the  diabatic  term,  note  that 
in  the  area  of  qualitative  verification  no  significant  dia- 
batic cooling  values  were  observed.   These  areas  are  limit- 
ed  in  extent  in  winter  and  spring  but  could  be  of  greater 
synoptic  consequence  in  the  summer  and  fall  seasons. 

The  effect  of  the  diabatic  term  as  it  appears  on  the 
forcing  function  side  of  equation  (21;)  gives  rise  to  figure 
(11),  applicable  at  approximately  850  mb,  but  actually  rep- 
resenting the  layer-mean  l3.   As  could  be  anticipated  from 
the  standpoint  of  heat  injection  into  moving  parcels,  the 
change  (inclusion  of  the  diabatic  term)  has  been  such  as 
to  expand  the  southern  rim  of  the  trough  which  extends  to- 
wards the  southwest  from  Iceland,   Similar  aspects  of  coher- 
ency between  the  resultant  co  of  this  study  and  the  trough 
movement  and  development  may  be  traced  out,  feature  by 
feature   as  one  proceeds  through  the  synoptic  sequences. 
For  example,  the  pronounced  trough  over  the  Atlantic  has, 
by  00Z  of  29  April,  become  oriented  North-South  just  east 
of  Greenland.   At  the  same  time  the  updraft  cell  has  taken 
on  this  same  orientation  just  west  of  the  same  trough,  so 
that  the  Li> -field  gives  some  confirmation  of  the  dynamic 
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processes  involved  in  these  map  changes. 

The  model  for  vertical  motion,  as  presented  by  this 
paper,  appears  to  be  quite  sucessful  and  future  plans  involve 
inclusion  of  a  more  realistic  frictional  term.  These  more 
realistically  computed  omegas  (extended  to  500  mb)  are  then 
to  be  used  for  feedback  to  yield  a  prognostic  z-field,  hour- 
by-hour.   The  ultimate  aim,  of  course,  Is  to  realize  a  small- 
er R.H.S.  error  as  the  motion  systems  progress  over  the  ocean, 
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